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52 I Abstract: A detailed model is constructed for the calculation of electron 

. ^ ' capture rates of some fp shell nuclei for situations prevailing in pre-supernova 

Vh ' and collapse phases of the evolution of the core of massive stars leading to 

supernova explosion. The model uses explicitly the Gamow- Teller strength 
function obtained through (n,p) reaction studies wherever available. The 
rates include contribution from the excited states of the mother as well as 
from the resonant states in equilibrium with the back reaction i.e. the beta 
decay of the daughter nucleus. Comparisons are made with the shell model 
results and the earlier calculations by Aufderheide et al. and Fuller, Fowler 
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and Newman. For the nuclei ^^Fe, ^^Mn and ^'^Ni with negative Q-values one 
observes large contribution from the excited states. 



Weak interactions play a significant role in the evolution of massive stars 
at the presupernova stage as well as at the subsequent core collapse phase of 
type II supernovae. The competition between electron capture and beta de- 
cay on fp shell nuclei (mostly with 50 < A < 70) at the late stages of silicon 
burning and the electron capture during the gravitational collapse, until neu- 
trino trapping sets in, determines the final lepton fraction inside the core, 
which in turn determines the shock energy at core bounce. Thus a care- 
ful evaluation of the electron capture and beta decay rates at presupernova 
densities and the electron capture rates at higher densities during the core 
collapse are of utmost importance for the supernova problem. The earlier 
detailed model of Fuller, Fowler and Newman [FFN,1] had all the important 
physics aspects incorporated in it. Later it was conjectured [2,3] that some 
nuclei with A > 60 not considered earlier may have significant impact on 
the deleptonisation and models for the calculations of the electron capture 
[4] and beta decay [5] rates for these nuclei were developed. Aufderheide et 
al [4] also identified sets of important nuclei for each density-temperature 
domain which need to be taken into account. Recently Ohu shell model cal- 
culations in the full fp shell with all the 4 orbits have been carried out for 
the evaluation of the electron capture and beta decay rates[6]. These rates 
are based on the allowed beta strength distributions using energy eigenvalues 
and wavefunctions obtained through very large matrix diagonalisation with 
a realistic fp shell interaction. 



On the other hand (p,n) reaction studies for a number of years and (n,p) 
reaction studies for the last decade or so carried out on a number of nuclei in 
this mass range have given us valuable information on the same allowed beta 
decay and electron capture matrix elements. It is always observed that both 
for beta decay and electron capture the experimental strength distributions 
are quenched compared to the best theoretical estimates. The models take 
account of this by including a quenching factor in the Gamow- Teller (GT) 
strength distribution, mostly a constant independent of energy. 

In an earlier work Kar, Ray and Sarkar [5] had calculated the beta decay 
rates for a number of important nuclei with A > 60 using a statistical model 
for the GT strength distribution including the effects of excited states and 
back resonances (coming from the equilibrium of beta decay and the inverse 
reaction of electron capture between the same two nuclei at typical high 
pre-supernova (preSN) temperatures). In this letter we extend that to the 
electron capture sector with the important improvement that the obseved 
(n,p) GT strength distribution is used directly for the evaluation of the rates 
for the ground state of the mother. So in this work we are concerned with 
the nuclei for which the (n,p) reaction results are available and we compare 
our results with those from the shell model calculations and other previous 
estimates. 

The electron capture rate on the nucleus {Z, N) is given by 



A ^^(p T Y) = In 2y^'^'^ l)exp[-Ei/kBT] fij{p, T, Fe, Qi, 
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G{Z,N,T) (ft) 
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Here p is the matter density, T the temperature, ks the Boltzmann con- 
stant and Ye is the electron fraction. G{Z,N,T) is the nuclear partition 
function at temperature T, Ei the excitation energy of the mother and Qij 
and {ft)ij are the Q- value and the ft- value between the i-th state of the 
mother and the j-th state of the daughter. 

Qij = Qg.s. + Ei — Ej (2) 

where Qg,s. is the ground state Q-value and Ei is the excitation energy 
of the i-th excited state of the mother and Ej is the corresponding energy of 
the j-th excited state of the daughter and 



ift), = , .:t;:^^. (3) 



(6250s) 

as for the ground state and the low-lying states there is no contribu- 
tion from the Fermi strength due to isospin selection rule. Bij{GT) = 
I Mij{GT^) I is the usual nuclear matrix element squared of the Gamow- 
Teller (GT+) operator. 

The phase space integral fij is of the form [7] 

1 f°° 
h = ^-^ / E{E^-miy/^iE-eofFe{Z, E)FD{E-i^e)FDGiE-fi,-eo)dE 

(4) 



with FD{x) = (1 + exp{x/kBT))^^ and FDC{x) = 1 — FD{x). rrie refers to 
electron mass, eo = —Qij and the lower limit of the integral Emin is the larger 
value between eo and nieC^- For the Coulomb correction factor Fc{Z,E) we 
use the form of ref [7] . 

For PreSN temperatures the reverse reaction on the daughter nucleus, i.e. 
beta decay of the daughter happens to be in equilibrium with the electron 
capture. For the ground state or a lowlying state of the daughter the state 
of the mother that is in equilibrium is not the ground state or other lowlying 
states, but one that has maximum overlap with the daughter state differing 
from it by a particle-hole excitation lying a few MeV higher in energy. We 
include the electron capture rate from this specific resonant state separately 
and write the electron capture rate as, 

>^fliiP. T, Fe) = \!^'{p, T, Fe) + X^siP, T, Fe) (5) 

where A^'"(p, T, Fe) is given by eq (1). A^£(p, T, Fe), the resonance con- 
tribution calculated using the model of Fuller, Fowler and Newman [1] is 
given by, 

A|£ = ln2 (6250s) -^(GVG"") x 

[(^)2 X exp{-Enes{GT)/KBT) x | Mgt \^ f{T, ^x,Qr,s{GT)) 
9v 

+ exp{-ER,s{F)/KBT) x | Mp ff{T,fi,QR,,{F)] (6) 



where C^ and G™ are the partition functions of the daughter and the 
mother, Ejies{GT / F) is the GT/F resonance energy in the mother and cal- 
culated using the model of FFN as described in ref [1] . 

QnesiGT / F) = Qg s.+EjiesiGT / F) for nuclei where the ground state sin- 
gle particle configurations of the daughter can get connected to the resonant 
state through GT/F transitions, otherwise QRes(GT/F) = Qg,s.- \ Mqt \ 
is the matrix element of Gamow- Teller operator for the reverse reaction i.e. 
beta decay [1] and is given by, 

I Mgt I' = ^„ E I M'a^ir, s) |' x ^^>\ (7) 

where n^ is the number of neutron particles in orbit 'r', n^ is the number of 
proton holes in orbit 's' (with degeneracy 2js + 1) and MQj^{r, s) is the single 
particle GT matrix element between orbit 'r' and 's'. Z„ is the quenching 
factor of the back reaction (taken as 0.6). | Mp \ = {N — Z) + 2 i.e. the 
difference in neutron and proton numbers of the daughter. ERes{F) is the 
excitation energy of the first higher isospin state (Tq + 1) in the mother 
{To = {N-Z)/2). 

For the Gamow- Teller strength distribution needed in the calculation of 
rates we use the results obtained from (n,p) experiments for the ground state 
of the mother. This is obtained by measuring the (n,p) reaction cross-section 
on the target and by performing multipole analyses on the resulting spec- 
tra. These experimental strength distributions are known to be substantially 
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quenched with respect to the theoretical predictions and also the shape of 
the energy distribution is different compared to the theoretical forms. The 
theoretical forms of strength distribution depend on the model space as well 
as the realistic interaction used. This is why the experimental strength dis- 
tribution is the best form to use for the rates wherever available even though 
there are sometimes large (about 20%) uncertainty in the extraction of the 
strength. For the nuclei ^^Fe [8], ^^Mn [8] and ^"^Ni [9] considered here we 
use directly the histograms of the GT strength per MeV from (n,p) studies. 
That is, we use the strength distribution as a function of energy instead of 
I Mij{GT^) I and replace the summation over the final states in eq(l) by 
integration over energy. 

For the excited states the rates are calculated using Gaussian forms for 
the GT+ strength distributions based on the spectral distribution theory 
arguments [10]. The same method was used in the beta decay sector in the 
model of Kar, Ray and Sarkar [5]. However as the excited states are within 
3-4 MeV of excitation we use the 'Brink hypothesis'and take for the centroid 
and width of the strength Gaussian the ground state values obtained from 
fits to (n,p) data. Thus the centroid of the strength Gaussian corresponding 
to transition from the i-th excited state is Ec{i) = Ec{g.s.) + Ei. The 
ground state strength centroid is fixed by the Sutaria-Ray expression [11] 
which gives a global fit to all known (n,p) data in fp shell as a function of 
A and A^. For the widths we use best fit values obtained by Sutaria-Ray for 



the specific examples [12]. 

Table 1 gives the calculated total rates for the nuclei ^^Fe, ^^Mn and ^''Ni 
with the four representative density grid points log pio =-2.5, -1.5, -0.5 and 
0.5 (pio is the density in lO^^g/cc) and five temperatures Tg (in 10^ K) = 
2,3,4,5 and 6. The rates within parenthesis are the ground state contribu- 
tions. The quenching factor in our theoretical model is taken as a parameter 
and fixed by the following procedure: For high density (log pio =0.5) and 
low temperatures (Tg =2) when the electron chemical potential is large to see 
the full GT strength distribution the quenching factor is adjusted to give the 
ground state rate identical to the one obtained by using the (n,p) data. This 
quenching factor is then used for all the excited states considered. The rates 
for all the three nuclei considered ( all with negative Q-value for e-capture) 
show a very rapid fall with decreasing density. This is because though the 
chemical potential (/ie) of the electrons is 11.98 MeV at log pio =0.5 (Tg=2,Fe 
=0.45) it becomes very small i.e. 0.804 MeV at log pio =-2.5 (Tg=2,y; =0.45) 
and at this lower density only a very small fraction of the electrons from the 
tail of the Fermi-Dirac distribution can overcome the Q-value and cause the 
capture. We also observe that for the highest density tabulated the contri- 
bution from the excited states is small (at most about 30% for the highest 
temperature). But at the lower densities when pe is too small to cause cap- 
ture on the ground state, the excited states still have a better chance for 
capture with Q-values higher (and sometimes even positive) and often give 



dominant contribution. On the other hand for the back resonances consid- 
ered though both Fermi and Gamow- Teller can in principle contribute , we 
find that their contributions even at log pio =0.5 is essentially zero (0.10 and 
0.00 for GT and Fermi respectively for ^"^Ni at Tg =6). For the nuclei ^^Fe, 
^^Mn and ^°Ni the GT back resonances are at excitations of 5.38, 6.71 and 
5.29 MeV respectively (according to the FFN prescription that we use) and 
these are too high to make any significant contribution to the rates. The 
back resonance excitation energies for the Fermi are all higher than 10 MeV 
and here again the exponential factor in eq (6) makes the contributions zero. 
In Table 2 we compare our rates with the rates calculated by other mod- 
els. We find that for the two even-even nuclei our rates are more than an 
order of magnitude higher than the ones by Aufdeheide et al. [4]. Actually 
Aufderheide et al. puts the total strength in a single excited state of the 
daughter nucleus with an effective 'log ft' and ref [6] observes that for the 
even-even nuclei it is put at too high an excitation. For the odd-A nucleus 
^^Mn however our values are within a factor of 2 higher compared to Aufder- 
heide et al.. Compared to FFN our values are within a factor of 2 on the 
higher side except for low density (p = 5.86 x lO^g/cc for ^°Ni) where the 
FFN value is an order of magnitude lower. But FFN used phenomenological 
log ft's concentrated in a giant resonance at a fixed energy whereas we use 
the observed (n,p) strength distribution. For ^^Fe at a density of lO^g/cc our 
ground state rate compares well with that of Martinez-Pinedo et a/. [6] but 
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our total rate is higher by a factor of seven. A more detailed comparison for 
other densities and other nuclei is needed to understand this. 

In conclusion we stress that for a proper comparison of the electron cap- 
ture rates by different methods a calculation which uses observed experimen- 
tal ground state GT strength distributions is important and is reported in 
this letter. We are in the process of extending our calculations to all the nu- 
clei for which (n,p) results are available. We also plan to look at the problem 
of identifying the sets of important nuclei for the preSN evolution for dif- 
ferent densities and temperatures using our rates and the nuclear statistical 
equilibrium. This may show important differences from the ones prescribed 
by Aufderheide et al. [4]. 
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TABLE 1 



e Capture Rates for Y^ = 0.45 
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Nucleus 


Temperature 
°ir 


logpio 


0.5 


-0.5 


-1.5 


-2.5 


56pg 


2x 10^ 
3x 109 
4x 109 
5x 109 
6x 109 


231.1 

(227.7) 

245.8 
(229.6) 

268.0 
(232.6) 

298.5 
(236.6) 

340.5 
(241.7) 


1.31 X 10-1 
(1.22 X 10-1) 

2.06 X 10-1 
(1.55 X 10-1) 

3.55 X 10-1 
(2.08 X 10-1) 

6.58 X 10-1 
(2.84 X 10-1) 

1.22 X 10° 
(3.19 X 10-1) 


1.65 X 10-^ 
(3.72 X 10-8) 

2.03 X 10-5 
(3.36 X 10-6) 

3.54 X 10-^ 
(4.52 X 10-5) 

2.41 X 10-3 
(2.63 X 10-^) 

9.45 X 10-3 
(9.73 X 10-^) 


8.21 X 10-11 
(1.13 X 10-11) 

1.14 X 10-^ 
(1.06 X 10-S) 

5.76 X 10-6 
(4.92 X 10-^) 

7.00 X 10-5 
(5.87 X 10-6) 

4.13 X 10-^ 
(3.57 X 10-5) 


55 Mn 


2x 109 
3x 109 
4x 109 
5x 109 
6x 109 


110.6 
(108.1) 

114.3 
(109.4) 

121.4 
(111.3) 

131.8 
(113.8) 

144.4 
(116.9) 


9.85 X 10-2 
(8.03 X 10-2) 

1.28 X 10-1 
(9.11 X 10-2) 

1.92 X 10-1 
(1.09 X 10-1) 

3.02 X 10-1 
(1.36 X 10-1) 

4.62 X 10-1 
(1.76 X 10-1) 


3.13 X 10-6 
(1.61 X 10-6) 

4.84 X 10-5 
(1.68 X 10-5) 

3.04 X 10-^ 
(7.76 X 10-5) 

1.15 X 10-3 
(2.48 X 10-^) 

3.17 X 10-3 
(6.57 X 10-^) 


1.46 X 10-9 
(4.00 X 10-1°) 

2.42 X 10-^ 
(5.41 X 10-8) 

4.40 X 10-6 
(8.54 X 10-^) 

3.05 X 10-5 
(5.56 X 10-6) 

1.31 X 10-^ 
(2.51 X 10-5) 
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TABLE 1 (contd.) 



e Capture Rates for Y^ = 0.45 



Nucleus 


Temperature 


logpio 


0.5 


-0.5 


-1.5 


-2.5 




2x lO'^ 


416.1 


6.66 X 10-1 


7.84 X 10-6 


4.99 X 10-9 






(415.5) 


(6.61 X 10-1) 


(3.66 X 10-6) 


(9.64 X 10-1°) 




3x 109 


426.7 


8.27 X 10-1 


2.51 X 10-^ 


1.66 X 10-6 






(418.7) 


(7.56 X 10-1) 


(6.84 X 10-5) 


(2.19 X 10-^) 


60Ni 


4x 10^ 


457.3 


1.25 X 10° 


2.45 X 10-3 


4.40 X 10-5 






(423.3) 


(8.96 X 10-1) 


(4.34 X 10-^) 


(4.74 X 10-6) 




5x 10^ 


514.8 


2.14 X 10° 


1.18 X 10-2 


3.65 X 10"^ 






(429.4) 


(1.09 X 10°) 


(1.64 X 10-3) 


(3.67 X 10-5) 




6x 10^ 


593.2 


3.45 X 10° 


3.55 X 10-2 


1.62 X 10-2 






(437.1) 


(1.35 X 10°) 


(4.58 X 10-3) 


(1.69 X 10-^) 



15 



TABLE 2 

Comparison of our electron capture rates with other model 

calculations 



Nucleus 


Density 
(x 10^ g/cc) 


Temperature 
(x 10^ ''K) 


Ye 


Our 

Rate 




Rates (s ^ 


)of 


Aufderheide 


FFN 


Martinez-Pinedo 












et al.[A\ 


[1] 


et al.[&\ 




5.86 


3.40 


0.47 


2.07E-06 


6.97E-08 


- 


- 


56pg 


10.7 


3.65 


0.455 


1.42E-05 


4.68E-07 


l.OE-05 


2.1E-06 




14.5 


3.80 


0.45 


3.97E-05 


1.31E-06 


2.81E-05 


- 




5.86 


3.40 


0.47 


2.17E-05 


1.49E-06 


1.09E-03 


- 


60Ni 


10.7 


3.65 


0.455 


1.25E-04 


7.64E-06 


- 


- 




14.5 


3.80 


0.45 


3.18E-04 


2.74E-05 


1.39E-04 


- 




33.0 


4.24 


0.44 


3.93E-03 


3.34E-04 


- 


- 




10.7 


3.65 


0.455 


1.52E-05 


9.23E-06 


- 


- 


55Mn 


14.5 


3.80 


0.45 


3.79E-05 


2.25E-05 


2.03E-05 


- 




33.0 


4.24 


0.44 


4.55E-04 


2.64E-04 


- 


- 
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